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Motivated b y new sounding-rocket wide- field polarimetric images of the Large Magellanic Cloud (Cole 
1999a), we have used a three-dimensional Monte Carlo radiation transfer code to investigate 



the escape of near-ultraviolet photons from young stellar associations embedded within a disk of dusty 
material (i.e., a galaxy). As photons propagate through the disk, they may be scattered or absorbed by 
dust. Scattered photons are polarized and tracked until they escape to be observed; absorbed photons 
heat the dust, which radiates isotropically in the far-infrared, where the galaxy is optically thin. The 
code produces four output images: near-UV and far-IR flux, and near-UV images in the linear Stokes 
parameters Q and U. From these images we construct simulated UV polarization maps of the LMC. We 
use these maps to place constraints on the star -f dust geometry of the LMC and the optical properties 
of its dust grains. By tuning the model input parameters to produce maps that match the observed 
polarization maps, we derive information about the inclination of the LMC disk to the plane of the sky, 
and about the scattering phase function g. We compute a grid of models with i = 28°, 36°, and 45°, 
and g = 0.64, 0.70^ 0.77, 0.83, and 0.90. The model which best reproduces the observed polarization 
maps has i — 36°l5 and g « 0.7. Because of the low signal-to-noise in the data, we cannot place firm 
constraints on the value of g. The highly inclined models do not match the observed centro-symmetric 
polarization patterns around bright OB associations, or the distribution of polarization values. Our 
models approximately reproduce the observed ultraviolet photopolarimetry of the western side of the 
LMC; however, the output images depend on many input parameters and are nonunique. We discuss 
some of the limitations of the models and outline future steps to be taken; our models make some 
predictions regarding the polarization properties of diffuse light across the rest of the LMC. 

Subject headings: polarization — methods: numerical — galaxies: individual (LMC) — ISM: dust, 
structure 



1. INTRODUCTION 



Polarimetric imaging provides a unique window on 



the 3 - dimcn3ional atructurc of aatrophyaical objects, and 



uum ultraviolet. This rocket-borne instrument has been 
flown three times to date, providing high-quality i mages 
of reflection nebulosity in the Pleiades open clust er ( |Gib- 



^ I thoroforo on the phyeical procoeeoe operating in a wide 



son, Holdaway, & Nordsieck 1995, Gibson 1997| , Gibson 



range of stellar and interstellar environments. The most 
important processes giving rise to interstellar polarization 
are scattering by dust, and transmission through aligned 
dust grains. 

Imaging polarimctry has been applied to many targets 



primary 



jf dotormining tho ooattoring propc 



fc Nordsieck 1999, in preparation), and ph otopolarimetry 
of Comet Hale-Bopp (Harris et al. 1997). Additionally, 
WISP obtained polarimetric images of the Large Magel- 
lanic Cloud (LMC); these observations represent t he first 
wide-field, UV p olarization images ever obtai ned (Nord- 
sieck et al. 199q ). Analysis of the LMC data (Cole et al. 



ti e s of dust, and to obtain geometric information on ex 
tended and complex sources (e.g., reflection nebulae, ac- 
tive galactic nuclei, and comets). 

The vacuum ultraviolet is an especially favorable wave- 
length regime for these studies; polarimetri c efficiencies are 



1999a, hereafter Paper I) found that the diffuse UV light is 



high, and polarized backgrounds are low (Nordsieck et al 



1993). Moreover, a relatively small number of bright stars 



emit the majority of VUV photons, greatly simplifying 
the accurate tracing of source-scatterer-detector geometry 
over the case in the optical and near infrared. 

The University of Wisconsin's Wide-Field Imaging Sur- 
vey Polarimeter (WISP) was developed to obtain the first 
wide-field astronomical polarization images in the vac- 



polarized at a 5-10% level, consistent with starlight scat- 
tered by dust; that the strongest source of illumination in 
the WISP field is the H II complex Nil; and that the UV 
starlight must account for most of the heating of diffuse 
dust in the LMC. 

In this paper, we report on our program to model the 
WISP polarization maps of the LMC using a Monte Carlo 
radiation transfer code to constrain the optical properties 
and scattering geometry of the dust in the LMC's diffuse 
interstellar medium. §1.1 and §1.2, respectively, describe 
the observational results which motivate this work, and 
the interpretive issues addressed by our models. 

In §2, we describe in detail the basic astrophysical ingre- 
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dients of our models, which are the distributions in space, 
size, and luminosity of the illuminating OB associations 
and scattering dust medium. In §3 we discuss the scatter- 
ing and polarizing properties of the dust grains and the 
way in which we parameterize the total amount of dust 
present. §4 describes the Monte Carlo engine of our radia- 
tion transfer code, which is innovative in its ability to track 
the processing of near-UV photons into thermal IR radia- 
tion by heated interstellar dust grains. In §5 the modelling 
algorithm is explained, relating the procedure by which 
we explored parameter space for the "best" models. §6 
presents the model images and polarization maps and the 
conclusions we can draw about the inclination of the LMC 
disk and the scattering asymmetry of its dust grains. We 
are careful to note the many shortcomings of this simple 
model, which nonetheless reproduces many of the observed 
near-ultraviolet and far-infrared properties of the LMC for 
a reasonable set of inputs. 



1.1. Observations: 



the Wide-Field Imaging Survey 
Polarimeter 



A 1?5 X 4? 8 area of the western side of the LMC was 
observed with the rocket-borne Wide-Field Imaging Sur- 
vey Polarimeter (WISP) on November 20, 1995. 4 x 80 
second exposures in an intermediate-band, near ultravio- 
let filter (A = 2150 A, AA = 300 A) were used to create 
intensity and polarization maps of the field. The observa- 
tions were centered at a = 04^ 59", 5 -67° 53' (J2000.0) 
and aligned roughly north-south. The WIS P instrument is 
described in detail in Nordsicck et al. 1993| ; the reduction, 
calibration, and analysis of the LMC flight data are given 
in Paper I. 

The minimum diffuse UV surface brightness, 5.6 ±3.1 
X 10~® erg s"""^ cm~^ A~^ Sr"-'^, is larger than any known 
stray light background, and is clearly due to light originat- 
ing within the LMC. The surface brightness of this diffuse 
UV background is correlated with areas of high H I column 
density and is linearly polarized at the ~10% level. This 
suggests that reflected OB starlight contributes at least 
half of the LMC's diffuse UV background. The ISM of the 
Large Magellanic Cloud apparently acts as a kiloparsec- 
scale reflection nebula in the near ultraviolet. Paper I 
found evidence for weak centro-symmetric scattering ha- 
los around some of the large OB complexes in the WISP 
field. The B2 complex (Martin et al. 1976), however, 
lacked such a halo; this was interpreted to mean that B2 
is located either within an H I hole or well above the plane 
of the LMC disk. 

1.2. Modelling Goals 

It is desirable to test the interpretation of Paper I; to this 
end we have undertaken to model the radiation transfer of 
ultraviolet photons from their origins in hot stars, through 
the dusty ISM of the LMC's disk, to Earth. Using these 
models we hope to determine whether or not the observed 
level of polarization is consistent with the reflection neb- 
ula interpretation. We also wish to determine the expected 
polarization pattern around B2 for a location within the 
disk; perhaps a non-detection of centro-symmetry is to be 
expected for this region. 

Using a specialized Monte Carlo radiation transfer code, 
we set out to determine whether or not reasonable values 
for the dust optical depth, scattering geometry, and dust 



grain optical properties can account for the WISP obser- 
vations. Under the assumption that the reflection neb- 
ula interpretation is correct, we can use the polarization 
properties of the model to place constraints on the dust 
properties and inclination of the disk of the LMC. 

2. MODEL INGREDIENTS AND ASSUMPTIONS 

2.1. Illuminating Sources 

The perfect model of the Large Magellanic Cloud would 
incorporate the luminosity contributions of every field star 
and star cluster into its input parameters. This is obvi- 
ously impractical, and so we must find a more tractable 
subset of objects with which to illuminate the LMC's 
dust. The star-formation rate of LMC field st ars has been 
roughly constant for the pas t «l-2 Gyr (e.g., Gallagher ei 



al. 1998; Westcrlund 1997 and references therein) . This 



recent activity has been accompanied by the formation 
of a large number of "blue populous" star clusters and 
OB associations; the young clusters of the LMC are both 
more frequent per unit field star mass and ind ividually 
larger than their Milky Way c ounterparts (e.g., Elson & 



Fall 1985; van den Bcrgh 1984) 



Data from the UIT instrument suggests that «75% of 
the flux from the LMC at A = 1500 A originates from stel- 
lar associations within the regions of nebulosity catalogue d 
by [Davies, EUiot fc Meaburn 1976| ([Parker et al. 1998|). 



For the western side of the LMC (observed by WISP), this 
interpretation holds true at 2150 A. In the WISP image, 
most of the wcll-detccted sources can be identified with 



OB associations (Lucke & Hodge 1970), or open clusters 
younger than «200 Myr. Clusters older than this, e.g., the 
massive 1 Gyr-old young globular cluster NGC 1783, are 
undetected in our image. Indi vidual supergiants among 
the field stars (Sanduleak 196£) can be detected, but are 
minor contributors to the total observed flux. 

We therefore chose to take as our illuminators the 122 



OB associations of Lucke & Hodge 1970, because a ho- 
mogeneous dataset of ultraviolet pho tometry at two wave 
lengths exists for the entire sample ( ^mith et al. (1987) 



hereinafter SCH). Due to the lack of a uniform sample of 
ultraviolet photometry, we have ignored the young open 
clusters in this first model; some of these clusters, e.g., 
NGC 1818, NGC 1755, and NGC 1711, contribute signifi- 
cant UV flux to the WISP image 



The posit ions of the OB associations from Lucke & 



Hodge 1970 were transformed onto the model's rectilin- 



ear coordinate system at a scale of 15' per grid unit. The 
scale was chosen in order to accomodate output images of 
the entire LMC, and the grid spacing is well-matched to 
the final, binned resolution of the WISP observations. 

The origin of th e coordinate system was chosen following 
Westcrlund 1990| to lie at 05^ 24™, -69° 50' (B1950.0); this 
corresponds to the ccntroid of o ptical light in the galaxy 
( deVaucouleurs & Freeman 1973). The distribution of OB 
associations in the LMC, as in the Galaxy, can be assigned 
some finite scale height above and below the galactic mid- 
plane. However, because this scale height is likely to b e 
smaller than the scale height of dust (Harris et al. (1997)), 
and we have no a priori knowledge of the relative positions 
of each association along the line of sight, we have forced 
the illuminating sources in our models to lie in the plane 
of the LMC disk. OB associations are not point sources, 
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having radii of w 15-150 pc ( |Lucke fc Hodge 197C ); for 
simplicity, we have modelled them as spheres. The radii 
of our iUuminators do not directly correspond to the op- 
tically defined dimensions of the Lucke & Hodge associ- 
ations, but were determined from the vacuum ultraviolet 
images of SCH. 

SCH photometered the entire Lucke & Hodge catalog 
at 1500 and 1900 A, using rectangular apertures that were 
matched to each association by hand; their Table 1 gives 
the total area of each of their apertures. We assigned radii 
to our spherical sources by setting their projected surface 
areas equal to the areas given by SCH in their Table 1. 
In just two cases (LH 15 and LH 77) did we deem the 
deviations from sphericality strong enough to warrant a 
more complex procedure. Both associations lie within su- 
pergiant shells on the northern side of the LMC disk. LH 
15, within LMC-1, is contained within the field of view 
of the WISP CCD image; LH 77, at the center of LMC- 
4, is quite bright and resembles a quadrant of a circle's 
circumference. These two associations were broken up ar- 
bitrarily into four identical sub-associations, which more 
closel y reproduced the visual appearance of these sources. 



directly. Harris et al. (1997) measured two reddening- 



Ncc r-UV luminosities were assigned to the sources based 
on the photometry of SCH at 1500 A (tois) and 1900 A 
(mig). We dereddeded the SCH photometry and apphed 
a correction for the difference in bandpass between their 
filters and t he WISP filt er at 2150 A. The reddening val- 
ues given by Lucke 1974 were broken down into LMC and 
foreground Galactic components; following SCH, the max- 
imum value of foreground reddening was taken to be E|;™^ 
— 0.07 mag. Any additional reddening towards the indi- 
vidual associations was attributed to dust within the LMC. 
Reddening values for each source were derived following 
the procedure outlined in Paper I, as were corrections for 
the differing bandpasses used by SCH and in Paper I. 

The derived extinction values wer e found to be in go od 
agreement with those published by Smith et al. 1990| in 
an erratum to SCH. The bandpass corrections ranged from 
—0.5 mag to -1-0.6 mag for the 122 Lucke-Hodge OB as- 
sociations. The corrected magnitudes were converted into 
monochromatic fluxes for the Monte Carlo photon gener- 
ator using the standard relation Fjjv = 10~°-^('"o+^^-^) 
(SCH). The source positions, radii, and luminosities are 
given in an appendix to this paper, in Table Al. 

2.2. Dust Distribution 

The LMC is a disk galaxy and we have chosen to rep- 
resent its dust density distribution using an exponential 
decay with radius a nd a hyperbolic secant law in height 
above the midplane (Binney & Tremaine 1987): 



p{r,z) = A exp( — ) sech^(- 



(1) 



where the constant A is set by the optical depth of the 
model (see §3 below), and the dust scale length and 
scale height are taken from the literature. Observa- 
tional estimates of the LMC's dust scale length were un- 
available, and so we set = 2.6 kpc (12 grid units ), the 



scale length of the old stellar population (Kinman et al 
19911 )^ 

The dust scale height, Zd, must also be estimated in- 
^We have assumed a distance to the LMC of 50 kpc 



free photometric indices for 2069 O and B stars in a 2.9 
deg^ area centered ^ 2?6 northwest of the optical center of 
the LMC. Using the distribution of reddening values they 
found that the data could be well-matched by a vertical 
distribution in which the dust has a scale height equal to 
twice that of the OB stars. Assuming the OB stars to lie 
in an extremely flattened disk, with scale height « 100 pc 
(c.f., [Oestreicher fc Schmidt-Kalcr 1995 for Galactic OB 
stars), we choose a dust scale height of w 200 pc, or 0.96 
grid units. 

As a flrst-order deviation from the smooth, azimuthally 
symmetric model dusk disk, we placed nine low-density 
cavities i nto the model, corresponding to the supergiant 
shells of Meaburn 198C. These shells were identified by 



the enhancements of Ha emission around their perime- 
ters and are also characterized by extremely low HI col- 
umn densities. They are thought to be roughly cylindrical, 
and "open-topped" (Westerlund 1997), but in our models 
they are defined by simple spherical cavities of low opti- 
cal depth. The cavities are placed in the midpl ane of our 
mode l galaxy, using the positions and sizes from 



Meaburn 



1980. We assign a density to the cavities by defining the 
near-UV optical depth Tc across the diameter of a cavity. 
Tc was derived from photometric measures of the redden- 
ing, FiB-v, of the OB associations lying within the bound- 
aries of the supe rgiant shells . These lie in the range 0.00 
<Eb_v'^0.12 ( [Lucke 197| ), less 0.07 mag of foreground 
reddening; we also assumed that roughly half of the ob- 
served reddening towards the OB associations was due to 
material in the near neighborhood of the stars and hence 
not a contributor to the optical depth of the cavity as a 
whole. We adopted a "typical" Eb-v of 0.01 mag, and 
assumed the OB associations to lie at the center of the 
spherical cavities; applying an LMC extinction law for the 
model's 2150 A photons, we set Tc = 0.1. The catalog of 
supergiant shell parameters is given in the Appendix, in 
Table A2; the cavities and illuminators are mapped out in 
Figure |l|. 

The inclination of the LMC disk to the plane of the sky 
remains a matter of some debate ( [Westerlund 1997| ). It 
has become clear that the east (30 Dor) side of the LMC 
is closer tha n the western (WI SP field) side. As shown in 
Table 3.5 of Westerlund 1997, both the inclination i and 



position angle line of nodes O are known only to a preci- 
sion of a few tens of degrees. Measurements of 9 scatter 
around a north-south line, and so we adopt 8 = 180° 
for simplicity. Geometrical methods applied to young and 
old stellar populations as well as neutral and ionized gas 
have yielded results varying between 25° <^ i <^ 48°. The 
expected magnitude and spatial variation of polarization 
depend strongly on the scattering geometry in the disk of 
the model galaxy, and so we consider three values for i 
in our models: 28°, 36°, and 45°. The WISP held, along 
the west side of the LMC, is tilted away from the Earth: 
as the inclination increases, photons must traverse larger 
path lengths through the absorbing dust layer in order to 
escape and be seen. 



3. DUST PROPERTIES 
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ninators (open circles) and cavities (dotted-line circles) in our model LMC. The positions nf ill nminatnrs 



distrihiition of illn 



are taken from Lucke & Hodge 1970, with radii and luminosities as described in §2.1. Cavity positions are taken from Meaburn 1980 Positions 
in equatorial coordmatcs nave been transformed onto a rectilinear grid with the origin at 5*^ 24™, —69° 50' (B195U.Uj. (Jttsets are given in 
degrees. North is up, East is to the left. OB Associations and supergiant shells in the WISP field (Paper I) are labelled with their Lucke-Hodge 
and Meaburn numbers, respectively. 
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Along with the scattering geometry, the optical proper- 
ties o f interstellar dust grains control the linear polariza- 



rect this val ue using the observed ultr aviolet (Pap er I) to 



tion cf stellar photons. Our models treat the scattering et al. (1997) found a mean B-band optical depth through 



proce ss using a standard Henycy -Greenstein phase func- 
tion (Henyey & Greenstein 1941), which depends on the 
albedo a, and asymmetry parameter g. The parameter 
g defines the probability for an incident photon to scat- 

ter through an angle 9: P{9) oc — — 



far- infrared ( DeGioia-Eastwood 1992) flux ratios. Harris 



"37 

(l+g2_2gcose) /2 

yields isotropic scattering, while g = 1 gives pure for- 
ward scattering; g < corresponds to backscattering. To 
model the polarization, we follow White (1979) in the ap- 
proximation P{9) w p^ax if^^- Pmax is the maximum 
polarization attainable in a single scattering event, for a 
scattering angle of 90°; the polarization of the scattered 
photon decreases for smaller and larger scattering angles. 

The total amount of dust in the model is described by 
the single optical depth parameter Teq. r^q is simply the 
optical depth of the model galaxy to a photon as it travels 
from center to e dge through the midplanc. 

White (lyvy / tabulated the scattering properties of the 



9 



Mathise^aL (1977)1 (MRN) Mil ky Way dust mixture. De- 
viations from the White (1979) values are to be expected 
for the LM C, which is in general more met al-poor than 
the Galaxy ( |Pagel et al. 1978| ; [Dufour 198| ). The LMC 
extinction curve shows a less pronounced 2175 A bump 
and a steeper rise into the far-UV than does the Galactic 
curve, attribut able to variations in du st grain size s and/or 
compositions (Nandy et al. 1981). Pci (1992) recalcu- 
lated the albedo of Magellanic Cloud dust using an MRN 
grain-size distribution with the relative contributions from 
graphite and silicates scaled to match the observed mean 
extinction curves. We adopt the Pei (1992) value, a = 0.66 
for A = 2150 A, for all our models; this is ~25 % higher 



than the o bserved alb edo of Milky Way dust (e.g. 
al. 1992 ). Pei (1992) did not include calculations for g or 



Witt e\ 



Pmax in his paper; we adopt the MRN-based value, Pmax 
= 0.31. Because the polarization of scattered starlight 
depends strongly on we compute families of models in 
which g is permitted to vary. 

In the ultraviolet, the phase function asymmetry pa- 
rameter g is poorly constrained by both models (0.1 <^ 
g ^0.7; e.g., Murthy & Henry 1994), and observations, 
(0.3 < g <;0.9; e.g., ^asseen &: Deharveng 1996| ), even for 
Milky Way dust. Variations in g lead to differences in 
the expected polarization patterns, allowing us to infer its 
value through comparisons of the model output to obser- 
vations. We chose to examine the moderately to strongly 
forward-throwing regime: g = 0.64, corresponding to the 
MRN value, and g = 0.90, suggested by observations o f 
reflection nebulae in the Pleaides cluster (Gibson 1997). 
Intermediate values were chos en at q = 0.70, 0.77, 0.83, 
close to the values derived by Witt et al. 1992 for the re- 
flection nebula NGC 7023. The most recent models of the 
Pleiades nebulosity (Gibson & Nordsieck 1999, in prepa- 
ration) also indicate a moderate value for g. 

For the LMC, the parameterization of dust mass us- 
ing Teq is problematic, because the LMC is observed 
nearly face on and thus the derived value of Teq is 
strongly geometry-dependent. We derive an ini tial opti- 
cal depth f rom observations of OB associations ( Harris e\ 
al. (1997)1 ), scaled by a geometric factor, and then cor- 



the face-on LMC tb — 0.98 x cos i. U sing the LMC mean 
extinction curve of Fitzpatrick 1986, the 2150 A optical 



depth to the LMC's midplanc is ttjv = ^ x 2.4tb = 1.2 



X cos i. The Harris et al. (1997) field is offset some 0.67 



radial scale lengths from the centroid of our dust distribu- 
tion, and so we adopt a pole-to-midplane central optical 
depth Tpoie = 2.4 X cos i. Teq is finally obtained by scaling 
Tpoie by the fiattening ratio of the disk. We adopt as our 
initial value: 



29.5 COS?. 



(2) 



This number was adjusted during the modelling proce- 
dure in order to match the observed ratio of WISP 2150 
A flux to IRAS 60 //m flux (see §5, below). 

4. RADIATION TRANSFER 

We construct model scattered light images with a Monte 
Carlo continuum radiation transfer code which accounts 
for multiple photon scattering, and predicts the spat ially 
resolved flux and polarization ( Wood fc Jones 1997 ). In 
the radiation transfer calculation, the dust plus gas mix- 
ture has albedo, a, and a scattering phase function approx- 



imated by the Heny ey-Greenstein phase function (Henyey 
fc Greenstein 194l| ) with asymmetry parameter, g (sec- 



tion 3). The code has been modifled from the axisym- 
metric models of Wood fc Jones to include a three dimen- 
sional distribution of illuminating sources (section 2.1) — 
this is crucial for modeling the UV scattered light pattern 
in the LMC. Additionally, we have added a new feature, 
where absorbed photons are not removed from the sim- 
ulation, but "re-radiated" isotropically from the point of 
absorption to form a "far-infrared" image. This is a first 
approximation for predicting the far IR emission from our 
simulations. We are in effect assuming that all the UV 
emission that is absorbed is reradiated at one wavelength 
where the dust is optically thin. In order to compare our 
models to the WISP and IRAS images of the LMC, we 
must include additional factors that account for the fact 
that our "reprocessing " technique does not enforce radia- 
tive equilibrium (section 5). 

5. MODELLING PROCEDURE 

With the dust distribution and OB association proper- 
ties held constant for all models, our goal is to find the 
combination of disk inclination and scattering phase func- 
tion that best reproduces the observed polarization maps 
of Paper I. The polarimetric images of the LMC are the 
result of the interplay of a large number of physical pro- 
cesses and properties; the vertical, radial, and azimuthal 
dust distributions, the dumpiness of the dust, grain albe- 
dos, the distances above or below the midplanc of each 
illuminating source, and the morphologies and densities of 
the supergiant shells all contribute to the observed polar- 
ization pattern. 

By tuning all of these quantities independently in our 
models, we could create a model whose output reproduces 
the observations precisely. However, such a model would 
almost certainly be non-unique, and its physical reality 
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would be highly questionable. We prefer instead to hold 
most of the model traits fixed and consider only the varia- 
tion in the parameters that are expected to most strongly 
influence the polarization maps. For example, one of the 
observables we insist the models match is the ratio of ul- 
traviolet to far-infrared light escaping from the galaxy. In 
reality, this obviously depends on the UV albedo a, but 
because tjr 1 and ruv > 1, it also depends on the opti- 
cal depth. We could trade off albedo and optical depth to 
tune this ratio precisely, but the model would suffer from 



severe non-uniqueness and could be driven to unphysical 



values of a, r, or both. Therefore we fix a (see Section 3), 
and adjust t until the model produces approximately the 
correct amount of infrared light. 

We consider a grid of 15 models (see Table 1), compris- 
ing 3 values of disk inclination i, and 5 values for scattering 
phase function g. For each model, we perform an initial 
run in order to determine the correction to our first guess 
at Teq. If the emitting dust is smoothly distributed and in 
radiative equilibrium, then the relation between absorbed 
starlight and thermal emission in the far-infrared depends 
simply on the dust optical depth. Hence we can attempt 
to match the amount of dust in our models to the true 
dust mass by a comparison of UV to IR fluxes; we correct 
our intial value of r^q by matching the model images in 
the UV and IR to the observed WISP 2150A and IRAS 60 
^m images. 

Rather than resort to large-scale averaging over ex- 
tended regions of the LMC, we choose instead to cali- 
brate our model optical depths using one well-measured 
region that lies within the WISP field of v iew (Paper I 



the N il complex (also known as DEM 3 4 [ Davies, Elliot 
& Me kburn 197611, Bl M^tm et al 1976( |, MClg |Mc(!!~ 
Brool s, fc: Batcnelor 1972[ |). Ni l lies at the sou thern edge 
of the supergiant shell LMC-1 ( Meaburn 1980| ), and con- 



tains a large H II region that is ionized by hot stars in OB 
associations LH9, LHIO, LH13, & LH14 ( |Lucke fc Hodje 
1970^ This permits us to make accurate comparisons be- 
tween our models and polarimetric observations. 

The correction is complicated because of the mismatch 
between models and reality. From the models, we compare 
monochromatic stellar photons at 2150 A to a monochro- 
matic far-infrared emission from dust which is heated by 
the starlight. From observations, we compare images in 
the A « 2150 A bandpass to emission at 60 ^m; real dust 
is heated by starlight of all wavelengths. Therefore the 
model IR/UV ratio cannot be immediately compared to 
the observations. In order to find the appropriate dust 
optical depth, we require: 



•^IRAS 
•^WISP 



UV UV IB. ^IR 
■ruv 



(3) 



where J-"iras and J^wisp are the observed fluxes in the 
IRAS 60 /im and WISP 2150 A bandpasses, and Tir and 
Tuv are the fluxes in the model output IR and UV images. 

r]^^ is a correction factor to account for the fact that 
our model dust grains are not in thermal equilibrium; 
they radiate their absorbed energy at a single, average, 
far-infrared wavelength whose flux equals the bolometric 



far-IR flux of the dust, rj^^ depends on the dust tempera- 
ture and the wavelength de pendence of the dust emis sivity; 
we adopt the values from DeGioia-Eastwood 1992 in her 
calculation of the ionizing flux in Nil. 

rf^^ is a similar correction that relates the dust heating 
by photons in the WISP bandpass to the total dust heating 
from light of all wavelengths. We calculate rf^ by inte- 
grating the light of Nil from 912 A to 3648 A using Ku- 
rucz model atr nosphcrcs, measurements of the in itial mass 
function from DeGioia-Eastwood 1992| (also see Parker e\ 
al. 1998), and weighting the spectral energy distribution 



by an LMC extinction law. 

^Yr corrects for the fact that our dust grains are not in 
radiative equilibrium; the model re-emits one far-infrared 
photon for each absorbed near-ultraviolet photon; to con- 
serve energy we must scale the model output images by 
the ratio of UV to IR mean photon energies, represented 

byes'. 

To tune the optical depth of our models, we checked each 
output model against equation ^; where the model ratio 
exceeded the the observations, we lowered Teq, and vice 
versa. By lowering the UV optical depth, we decreased 
the number of reflected photons that are subsequently ab- 
sorbed, and hence the relative amount of far infrared emis- 
sion. As a result of the tests, we adopted the parameters 
shown in Table 1 for each inclination; no systematic trend 
with g was apparent, although the values of showed 
a scatter of ±10% around the mean value at each inclina- 
tion. Due to difficulties with the IRAS zeropoint calibra- 
tion, and uncertainty in the WISP zeropoint, it is difficult 



to relate our model parameter 



to a physical fiux ratio. 



Because of the large number of photon sources, the num- 
ber of output images, and the large number of pixels in 
each output image, large numbers of photons were required 
in order to obtain significant signal-to-noise to measure the 
polarization of the diffuse UV light. For each of the 15 per- 
mutations of i and g, we computed a "low" signal-to-noise 
model with 10^ photons propagating through the model 
galaxy. These models provided sufficient information to 
identify the models that matched the data well enough to 
merit a more detailed look. 

We re-ran the best fitting model with 10^ photons in 
order to more accurately assess the mean level of polar- 
ization and track the variation in polarization level and 
position angle across the image. As we began to write up 
these results, we continued to let the model run in order 
to create the highest possible signal-to-noise in the output 
images. 



6. RESULTS & DISCUSSION 



6.1. 



Model Images & Polarization Maps 

Our models produced as their output a set of four im- 
ages for each run: the ultraviolet and infrared flux, and, in 
the ultraviolet, images of the linearly polarized flux Q and 
U. For analysis, the output data were converted into FITS 
format and examined within the IRAF^ suite of tasks. 

Model images of the WISP-observed portion of the LMC 
are shown in Figure |^. From left to right, we show the 
model UV image, the WISP 2150 A image, the IRAS 60 
image, and the model IR image. The far left and right 



^IRAF is distributed by the National Optical Astronomical Observatories, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Table 1 

Assumed, derived, and output parameters for our Monte Carlo models. 





(deg) 


'eg 




5yi'^ (arcmin) 


0.64 


28 


24.8 


12.2 ±0.3 


0.6 ±0.6 


0.64 


36 


22.4 


12.4 ±0.3 


10.2 ±0.6 


0.64 


45 


20.1 


12.8 ±0.3 


16.2 ±0.6 


0.70 


28 


25.9 


12.7 ±0.3 


-0.6 ±0.6 


0.70 


36 


23.3 


12.6 ±0.3 


9.6 ±0.6 


0.70 


45 


19.9 


12.9 ±0.3 


15.6 ±0.6 


0.77 


28 


25.5 


12.9 ±0.3 


0.0 ±0.6 


0.77 


36 


23.5 


12.8 ±0.3 


8.4 ±0.6 


0.77 


45 


19.8 


13.0 ±0.3 


16.2 ±0.6 


0.83 


28 


25.3 


13.1 ±0.3 


-0.6 ±0.6 


0.83 


36 


24.5 


13.4 ±0.3 


6.6 ±0.6 


0.83 


45 


20.3 


13.2 ±0.3 


15.6 ±0.6 


0.90 


28 


26.5 


13.4 ±0.3 


-1.2 ±0.6 


0.90 


36 


24.4 


13.7 ±0.3 


7.2 ±0.6 


0.90 


45 


20.8 


14.1 ±0.3 


13.8 ±0.6 


Observations' 




12.6 ±2.3 





"^Scattering asymmetry parameter, see §3. 

^Inclination angle of LMC disk, see §2.2. 

■^central optical depth in the plane. 

^mean percentage of linear polarization. 

"offset in polarization centers of symmetry, west of 
illuminating source. 

ffrom Cole et al. 1999a. 
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MODEL POLARIMETRY 



panels of Figure g show the co-added results of multiple 
Monte Carlo runs, amounting to a total photon count of 
4.1 X 10^. Only a portion of the image is shown, to fa- 
cilitate comparison to the WISP data (middle left panel 
of Figure ||) . The dynamic range of the UV model is sig- 
nificantly smaller than the 2150 A data, because of the 
uniformly high optical depth in our models which results 
in nearly constant attenuation across the field. Many of 
the features are well-reproduced by the models, including 
the prominent Nil and B2 complexes which particularly 
dominate the IRAS 60 fim image. The supergiant shells 
are less visible than expected in the model IR image; it 
lacks the "holes" visible in the 60 fim data. However, 
their effect can be noted by comparison of the UV and IR 
appearance of the LII15 association at the upper right of 
the panels in Figure ^: LII15 is far brighter in the UV 
than the IR, a consequence of its location within the shell 
LMC-1 (see Figure |l|). 

Despite the obvious, and expected, lack of fine-scale 
structure in our model, we nonetheless are able to approx- 
imate the large-scale luminosity distribution in these two 
wavelength regimes. The most obvious failure of our model 
to reproduce small-scale structure in the dust distribution 
occurs at the association LH12 (see Figure ||). LH12 is 
an intrinsically bright association (SCH) that is highly 
reddened (Lucke 1974) and so appears faint to observers. 
However, LH12 lies within the LMC-6 supergiant shell; 
in our simple model, this drastically reduces the amount 
of obscuring dust along the LH12 sightline, causing us to 
dramatically overestimate the observed brightness of the 
association. A similar effect is discernible for the more 
southerly association LH5, highly reddened despite its po- 
sition near the edge of the LMC-7 supergiant shell. 

We created images of the degree of linear polarization, 
p, and position angle 9: 



P = 



/(Ql 



U2) 



P 



0.5 



0.5 arctan 



(4) 
(5) 



We created polarization maps with vector length pro- 
portional to P and position angle equal to 9; our highest 
S /N map is shown in Figure ^. The smooth dust distribu- 
tion and finite number of illuminating sources account for 
the extreme regularity of the model polarization map. 

In order to constrain the disk inclination of the LMC 
and the phase function g of its dust, we compared the 
model polarization maps to those produced by Paper I 
from the WISP observations, shown side-by-side in Figure 
The WISP observations, binned into 6' pixels to in- 
crease signal-to-noise, showed some evidence for the pres- 
ence of centro-symmetry about the brightest OB associa- 
tions, but were hampered by the shortage of photons far 
from these bright regions; the effects of inhomogeneities in 
the scattering medium and a diffuse starlight component 
clearly dominate the appearance of the observed polariza- 
tion maps. 

Our points of comparison included the mean P and 9, 
the shape of the high-polarization tail of P, the degree of 
centro-symmetry around the brighest regions (Nil [Bl] in 
North and N79 [B2] in the South), and the offsets of the 



polarization symmetry centers from the central OB asso- 
ciations. We also looked for variations in mean P and 9 
across the 1?5 x 4?8 field of view, but found no statistically 
significant differences to distinguish the various models. A 
The comparison was not straightforward because of the 
systematic effects which dominate the errors in the ob- 
served polarization maps. In particular, the mean value 
of P was determined quite carefully: the WISP maps have 
been cleaned of marginal polarization detections where P 
< 2.5crp, and become incomplete for polarizations less than 
10%. Therefore we applied a linear incompletess correc- 
tion to the models, such that the model polarizations are 
weighted progressively less between 10% and 4%, and all 
polarizations smaller than 4% are ignored. This mimics 
the observational bias towards detection of regions of high 
polarization. 

The observed polarization maps, biased by incomplete- 
ness toward the detection of high polarization regions, 
showed a much higher scatter in polarization values than 
did the models. The observations also show a patchier 
distribution of P; pixels with polarizations in the 5-10% 
range are frequently juxtaposed with w20% polarized re- 
gions. This indicates that the effects of small, optically 
thick clumps in the ISM are strongly influencing the scat- 
tered light component of the diffuse UV background. How- 
ever, when we binned the models to the resolution of the 
observations and corrected for incompleteness at P < 10%, 
our models were consistent with the observed mean level 
of polarization across the WISP field. Scattered light from 
OB associations is indeed likely to account for a large frac- 
tion of the diffuse ultraviolet background in the LMC. 

6.2. Disk Inclination and the Dust Scattering Phase 
Function 

We found that our model polarization maps were sensi- 
tive to variations in g and i; see Figure |^. As g increased 
from 0.64 to 0.90, the mean polarization (P) increased 
from Ril2.5% to «13.7%. The disk inclination manifested 
itself most noticeably in the distribution of polarization 
vectors around bright sources, e.g.. Nil. For a face-on 
disk of scattering material, the polarization vectors form 
a centro-symmetric pattern about the illuminating source; 
the inclination introduces an asymmetry which shifts the 
center of the distribution away from the illuminator, per- 
pendicular to the disk's line of nodes. As the inclination 
was increased, the offset in symmetry center of the polar- 
ization patterns (hereafter referred to as Jx) varied from 
0'.36 eastward, at 28°, to 16f2 westward, at 45°. These re- 
sults are summarized in Table 1, and plotted in Figure ^ 
in which (P) and 6x are plotted for each model. We have 
also plotted (P) and 6x for the observed polarization map, 
with the associated error bars, ^x tends to zero at incli- 
nations of <,30° due to the competing effects of the disk 
inclination and the radial drop in dust density from east 
to west across the field; a plane parallel slab of scattering 
dust would show 6x = only for i — 0. 

Interpolating in inclination, and adopting the uncer- 
tainty in our observational determination of ^x from the 
binned WISP data, we find that the LMC's disk is inclined 
at 36° ^5 to the plane of the sky. Insofar as the western 
side of the LMC resembles our models, this is a direct, 
geometric determination of its inclination. Although in- 
consistent with recent kinematic determinations (e.g., Kim 
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Fig. 2.— Model output images for i = .'^fi" rj = n 70 and 4.1 X 10^ total photons. Prom left to right: model UV image; WISP 2150 A 
image (Paper I); IRAS 60 /jm image (Schwering 1989); model IR image. North is up, East is to the left. The differences between model and 
observations are attributable to the highly non-unitorm distribution of stars and dust in the LMC. 
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Scale: — 10% polarization 
Angular Scale: 

1 degree 

Fk;. 3. — Model UV image with the polarization vectors derived from the Stokes Q and U model images overplotted. Because the only 
illumination derives from the OB associations, and the dusty scattering medium is smoothly distributed, the polarization vectors present 
clear, regular centrosymmetric patterns about the illuminators. North is up, east is to the left. 
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Fig. 4. — left: model UV image with polarization vectors overplotted; right: WISP 2150 A image with polarization vectors. The polarization 
and angular scales are shown in the figure. 
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Fig. 5. — Model polarization maps for the area surrounding the H II complex N 11 demonstrate the dependence on inclination and phase 
function asymmetry. The actual 2150 A image of the region, 1° on each side, is shown at upper left. Polarization maps for the same area are 
shown, arrayed by inclination and g value. The smallest polarization vectors plotted are 1%, the largest 35%; observed values around N 11 
range from to ~10%. 
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et al. (1998), who found i = 22° ±6°), it is squarely in 
agreement with th e average of rehable values tabulated by 
Westerlund 1997. If the dust scale height is larger than 



we have assumed, then the inclination required to yield a 
given (5x is reduced; a significant change in Zd would be 
required due to the concentration of scattering dust close 
to the midplane of the disk. 

It is more difficult to draw conclusions regarding g. As 
seen clearly in Figure ^, the effect of high g is mainly to 
increase the mean level of polarization. This ran counter 
to our expectations, for a high value of g should produce 
smaller scattering angles, on average, resulting in polariza- 
tions much smaller than Pmax ■ The counterintuitive result 
is due to our simulation of the observational biases in the 
data of Paper I (see section 6.1); the high g models did ini- 
tially show lower levels of polarization, but were cleaned 
more severely by our rejection of low P pixels. 

The low signal-to-noise of the WISP data and the larger 
than expected polarization fluctuations result in uncom- 
fortably large errors in our determination of (P). At the 
3(7 level, we cannot rule out any value of g between and 1. 
At the Icr level, we find g = 0.7 ±0.3. The simple compari- 
son of (P)mod to (P)obs throws away information about the 
distribution of (P) within the images: we found that for g 
= 0.83-0.90, our models exhibited a strong high-P tail ex- 
tending to «30%. This tail was weak in the observations, 
and suggests that such high values of g are less likely than 
would be inferred from the formal error distribution. How- 
ever, as noted in Paper I, the highest polarization pixels 
typically show low flux values, and hence relatively large 
errors in both P and 6. Low signal-to- noise might therefore 
introduce a bias against the high-g models. 

Our determination of g is of comparable p recision to 
value s for t he di ffuse Milky Way dust, e.g. [Murthy "& 
Henr} ] 1994| and | asseen fc Deharveng 199q ; because of 
the complex scattering geometry it is far less precise than 
determinations ba sed on Galactic reflection nebulae, e.g., 
Witt et al. 1992. Our value depends critically on the 



smoothness of the dust distribution as well as the relative 
positio ns along the linc-of-sig ht of illuminators and scat- 
terers ( |Witt fc Gordon 1996i Gibson & Nordsicck 1999, 
in preparation). Also, as shown in Figure ^, the true dis- 
tribution of ultraviolet luminosity in the LMC contains a 
non-negligible component due to stars not in the Lucke- 
Hodge associations. This additional source of direct, un- 
polarized, light would dilute the scattered light and hence 
reduce (P) if included in the models. In such a model, 
higher values of g would be required to account for the 
observed level of polarization. 

6.3. Individual associations 

Paper I lists nine UV-bright regions which appear 
to support scattering halos with the expected centro- 
symmetric pattern; we explore here the degree to which 
the simple model is able to reproduce these features of the 
polarization map. 

6.3.1. NGC 1755, NGC 1711, N 186 

Two open clusters and a modest-sized star-forming re- 
gion, these objects show some evidence of scattering halos; 
however, they are not included in our models. 



6.3.2. LH 15, N 11, LH 12 

These bright OB associations support scattering halos. 
The models reproduce them quite strongly. The northern 
associations LH 15 and Nil dominate the local UV radia- 
tion field, and are well detected in the WISP observations. 
LH 12 is far brighter in the model than is observed, a con- 
sequence of the strong deviation from smoothness of the 
surrounding dust structures (see §6.1). The observed scat- 
tering halo around LH 12 is quite weak, indicative of the 
contribution of increasing field star density in the southern 
half of the WISP field (Paper I). 

6.3.3. LH I LH 25, LH 16-17-20 

These regions show very weak evidence of scattering ha- 
los in the WISP image. In the models, there is little ap- 
parent indication of centro-symmetry about these associ- 
ations. We note that each of these associations lies very 
close to the edge of the WISP image. Moreover, they lie 
near or within the scattering halos of brighter associations. 
These must be considered marginal detections. 

6.3.4. B2 

The B2 complex is made up of associations LH 1, 2, 5, 
and 8; its large OB star population led us to expect the 
presence of a strong centro-symmetric scattering halo. To 
the contrary, the paper I analysis of the WISP observations 
found no such halo. Possible explanations were suggested: 
the location of B2 above the plane of the dust layer, loca- 
tion of B2 within a large H I hole, or a possible bias against 
detection of scattering halos larger than w40'. Our model 
polarization maps show a centro-symmetric halo around 
B2; however, comparison to the comparable associations 
Nil and LH 15 showed it to be weaker than the halos of 
the northern associations. This is due, in the models, to 
the extnded (non-point-source) size of B2, as well as the 
contribution of light from LH 12 and other associations 
beyond the WISP field of view. In addition, the observa- 
tions are degraded by the presence of high field star den- 
sity in the southern WISP field, as well as the existence 
of a significant hole in the H I distribution (S. Kim, pri- 
vate communication). Part of this H I hole corresponds 
to the Ha supergiant shell LMC-7 (Meaburn 1980), which 
coincides with LH 8. We conclude that the Paper I non- 
detection of a scattering halo around B2 is astrophysical 
and not due to bias in their analysis, but the placement 
of the complex above the plane of the dust is not required 
by the observations. 

6.4. Limitations & Future Work 

This suite of models represents our first attempt to 
model the radiation transfer of polarized light through a 
galaxy from a large number of discrete sources within a 
non-uniform dust layer. These models are able to repro- 
duce the general morphology of ultraviolet and infrared 
images of the LMC; taken together with observed polariza- 
tion data, they yield astrophysically interesting constraints 
on its inclination, and show consistency with dust proper- 
ties expected from observations of the Milky Way. 

However, the models provide a greatly simplified pic- 
ture of the true structure of the LMC. Future adaptations 
of the radiation transfer code will address many of the 
simplifications; for others, additional observational mate- 
rial is required in order to refine our input parameters. A 
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Fig. 6. — Model images created with differing values of i and g separate themselves in a plot of mean percentage of polarized flux against the 
offset between the observed scattering halos and illuminating sources. Sx is shown in arcminutcs. The solid square shows the observational 
data with errorbars (Paper I). A typical crrorbar for the modcl-dcrivcd points arc shown in the lower right corner. Each nearly vertical column 
of symbols is labelled with the disk inclination; the legend relating symbol type to g is at lower left. 
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more advanced version of the Monte Carlo code (dubbed 
"galaxy on a grid" ) , will allow the specification of dust and 
luminosity density at each individual point in the model 
grid; this flexibility will permit the exploration of arbitrar- 
ily complex geometries in future work. 

The most severe drawback of our model is the adoption 
of a smooth dust distribution that lacks o ptically thick 



clum ps. As noted by e arlier authors (e.g., Witt, Thron- 



son ^ Capuano 1992| ; [Witt fc Gordon 1996| ), the pres- 
ence of small-scale, dense knots of absorbing material can 
greatly alter the emergent spectral-energy distribution of 
a galaxy's light. In this case, inferred values of g are incor- 
rect, and our ability to distinguish d ust optical pr operties 
from the scattering geometry is lost ( |Gibson 1997 , Gibson 
& Nordsieck 1999, in preparation). Since the distribution 
of small-scale dust knots in the LMC is unknown, we have 
little recourse but to adopt a relatively homogeneous dis- 
tribution of dust. The obvious effect of this approximation 
is that our models lack the large polarization fluctuations 
seen in the observations, presumably attributable to the 
presence of the patchy distribution of optically thick dust 
clouds. 

A related drawback is the unknown line-of-sight distri- 
bution of illuminators relative to the dust. As shown by 



Witt, Thronson & Capuano 1992, the relative amounts 



of scattering and absorption depend sensitively on where 
the illuminating sources lie relative to the dust. Along 
the same lines, we have modelled the supergiant shells as 
spheres, whereas evidence suggests they more closely re- 
sembl e cylinders. One obvious effect of this mismatch is 



the pi esence of a foreground "haze" of IR emission above 
the shells in our models, which is not present in the IRAS 
data (see Figure ||) . 

The dust scale length in our models is highly uncertain; 
if we had chosen to adopt the scale length of H I inferred 
from the maps of |Kim et al. (1998)| , our derived would 
have been some 40% smaller than the value we used. This 
in turn would have required a higher r^q in order to repro- 
duce the UV— IR color of the diffuse light at the position of 
the WISP observations. Future models, using the galaxy 
on a grid system, would be set up to more closel y match 
the c omple x column-density variations seen in the Kim e\ 
al. (1998) H I maps, and therefore obviate the necessity 
of parameterization using r^. 

Although our scattering medium is unrealistically 
smooth, our stellar sources suffer from the opposite prob- 
lem: we have considered only emission from large OB as- 
sociations, ignoring the contributions of young, massive, 
open clusters and the field star population. This yields a 
model with a very highly clumped luminosity distribution; 
because our models ignore the direct stellar contribution 
to the diffuse ultraviolet light, they will tend to produce 
higher levels of polarization that would otherwise be ob- 
served. For this reason, our predicted value of g is likely 
to be skewed towards lower values than a more complex 
model would produce. In future work, it will be desir- 
able to add a smoothly varying component to t he starlight 



taken for example fr om the large-scale maps of Mauchcrat- 
Joubelrt et al. 19801. 



6.5. Predictions 

Our models extend over the central «10° x 10° of the 
Large Magellanic Cloud, roughly 14 times the area covered 



by the WISP observations. This allows us to predict the 
general pattern of polarization that might be seen across 
the face of the LMC at near- ultraviolet wavelengths. Our 
models have been tuned to reproduce the optical depths 
and polarizations of the WISP field; a test of their validity 
and generality would be a comparison of our predictions to 
future ultraviolet polarimetric datasets across rest of the 
Large Cloud. 

In Figure |^, we show the ultraviolet polarization vectors 
from our best model for the entire LMC. An optical image 
of the galaxy (Sandage 1961) is plotted to provide orienta- 
tion and scale. If our model is valid. Figure |^ should pre- 
dict the pattern of polarization across the face of the LMC. 
The predictions of the model are most uncertain in the re- 
gion of the bar, which contains a high-surface brightness, 
intermediate-age population of field stars and numerous 
youn g open clusters not acou nted for in our models (see, 
e.g., Hodge & Wright 1967). We comment upon some 
regions of interest that may be likely future targets for 
polarimetric study. 

We see that the 30 Dor region itself shows low levels 
of polarization, although it is responsible for much of the 
scattered light within a kiloparsec or more. The bar is 
expected to show low levels of polarization (^ 5%), with 
most of that due to 30 Dor at the eastern end. The bright 
H II region N51 (north of the central bar) should produce 
a centrosymmetric pattern similar to that of Nil (in the 
WISP field), although it will be weakened on the eastern 
and southern sides by radiation from 30 Dor. Finally, the 
southe rn spiral arm seen in the H I maps of ^ im et 



(1998), beyond the southern limit of Figure t], is predicted 



to scatter light from as far north as 30 Doradus, giving rise 
to faint diffuse light that is polarized at the 20-30% level. 
It is expected that these models overpredict the mean lev- 
els of polarization because of the neglected contribution of 
direct light from field stars. 

7. SUMMARY 

Using a Monte Carlo radiation transfer code, we have 
modelled the observed ultraviolet polarization maps of the 
LMC obtained with the WISP instrument. Our code fol- 
lows the tracks of stellar photons from their origins within 
OB associations, through a smoothly distributed exponen- 
tial dust disk containing low-density cavities. By account- 
ing for far-infrared thermal emission from heated dust, we 
are able to parameterize the total amount of dust present 
by its optical depth. 

Dust-scattered starlight gives rise to linear polariza- 
tions; the magnitude and position angle of the polariza- 
tion vectors allow us to derive information regarding the 
scattering geometry of stars-|-dust within the LMC and 
the optical properties of Magellanic Cloud dust. We con- 
sider three disk inclinations between 28° and 45°, and five 
values for the phase function asymmetry parameter g, be- 
tween 0.64 and 0.90. We derive: 

• The inclination of the disk of the LMC to the plane of 
the sky is 36° This is in agreement with other results 



( Westerlund 1997), but not w ith a recent kinematic deter- 
mination (Kim et al. (1998)). Our determination contains 
a dependence on the dust scale height; we have assumed 
Zd — 200 pc, but smaller values would imply larger incli- 
nations, and vice versa. 

• The most likely value for g of 0.70; the uncertainty in 
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Scale: — 1 09b polarization 



Angular Scale: 



I degree 



Fig. 7. — Polaximetric predictions for the LMC: our model UV polarization vectors are overlain on an optical image (A Ri6600 A) of the 
galaxy from Sandage (1961). The WISP field of view is at the right. North is up, East is to the left. 
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the observations does not permit us to rule out any value 
for g. Higher signal-to- noise data would sharpen our es- 
timate, but a more precise method for estimating g from 
our models is needed as well. For g above 0.77, the models 
produce a higher fraction of highly polarized {p ^ 20%) 
pixels than are observed. However, the neglect of direct 
light from field stars probably leads us to understimate g. 
Our value g k, 0.7 is consistent with values derived from 
ultraviolet surface photometry of t he Galactic reflection 
nebula NGC 7023 ( |witt et al. 1992| ). 

• Our best model predicts that scattered light from 30 
Doradus dominates the eastern side of the LMC's diffuse 
UV radiation field; the presence of this "mini-starburst" 
region may be felt as far south as —72°. Data for this 
region would be of great value in constraining the LMC's 
structure and geometry. 

Our models are the first attempt to replicate the ob- 
served images and polarization maps of the Large Mag- 
ellanic Cloud. It is encouraging that our model is able 
to approximately match the observations of Paper I for 
a reasonable set of input parameters; we do not expect 
that we have found a unique solution for i, g, or the geo- 



metric parameters that go into the model. The modelling 
procedure should be applicable to more general problems 
in the interpretation of galaxy polarization maps. How- 
ever, a large number of simplifying approximations and 
assumptions have been made. The most serious of these 
is the smoothness of the adopted model dust distribution. 
Future work will be able to take into account the complex- 
ity of the LMC's H I distribution, which is dominated by 
flocculent spiral arms, super giant shells, and small scale 
filaments (|Kim et al. (1998)|). 
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APPENDIX 



MODEL GEOMETRY 



Large OB associatons are assumed to dominate the diffuse ultraviolet radiation field of the model galaxy. The source 
list has been taken from the catalog of Lucke & Hodge (1970), and associations are identified by their LH number. Their 
radii and fluxes have been derived for our model purposes in §2.1. Our adopted source parameters have been listed in 
Table Al. Source positions are given in terms of the rectilinear model coordinate grid, where x and y are the east-west 
and north-south offsets, respectively, from the optical center of the LMC bar at 05*^ 24™, —69° 50' (B1950.0). The source 
radii do not replicate the true physical sizes of the LH associations, but are the effective radii of circular regions of the 
equivalent areas of the associations (see §2.1). 21 of the associations lie within the area observed by WISP (Paper I). 

The smoothly distributed exponential disk of dust in our models has been seeded with cavities that approximate the 



extremely low density supergiant shells investigated by Meaburn 198C. The supergiant shells were identified by Ha 



emission from their limbs, and roughly correspond to areas of low HI column density (Kim et al. 1998). While the real 
supergiant shells are approximately cylindrical, our code as currently implemented allows only spherical cavities. The 
resulting shells are overlain by the high-z tail of the vertical dust distribution an d thus fail to p recisely reproduce the 
infrared morphology of the shells. The cavity parameters, with identifications from Meaburn 1980|, are given in Table A2. 



In each case, the optical depth at 2150 A across the diameter of a cavity has been set to 0.1 (see §2.2). Three of the shells 
are contained within the observed WISP field (Paper I). 
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Table A1 

Source parameters for Monte Carlo models. 



LH^ 




x'^ 






flux'* 


LH 


X 


y 




flux 


LH 


X 


J 


Tcld 


flux 


1* 


2 


312 


-0 


250 


4.05 


4.91 


42 


-0 


050 


-0 


300 


6.00 


14.45 


83 


-1 


580 


3.240 


3.45 


5.01 


2* 


2 


250 


-0 


200 


2.10 


3.40 


43 


-0 


100 


3 


750 


4.35 


5.75 


84 


-1 


532 


2.380 


4.65 


26.46 


3* 


2 


467 


1 


732 


4.35 


0.87 


44 


-0 


200 





267 


6.00 


0.42 


85 


-1 


500 


0.500 


2.70 


14.20 


4* 


2 


467 


2 


132 


4.35 


6.26 


45 


-0 


200 


3 


480 


8.40 


12.91 


86 


-1 


500 


1.820 


2.55 


0.79 


5* 


2 


100 





000 


4.05 


8.69 


46 


-0 


250 


-0 


080 


3.30 


2.09 


87 


-1 


450 


-0.240 


4.95 


46.07 


6* 


2 


267 


1 


932 


4.35 


3.98 


47 


-0 


267 


1 


380 


4.50 


92.73 


88 


-1 


532 


1.732 


2.10 


0.46 


7 


1 


700 


-2 


160 


2.85 


0.08 


48 


-0 


280 


1 


425 


2.85 


7.37 


89 


-1 


500 


0.410 


5.85 


93.84 


8* 


1 


920 


-0 


250 


6.75 


10.39 


49 


-0 


350 


1 


250 


2.85 


4.87 


90 


-1 


500 


0.200 


2.85 


19.17 


9* 


2 


240 


2 


600 


3.00 


34.88 


50 


-0 


450 


-1 


932 


6.75 


1.10 


91 


-1 


650 


2.780 


2.55 


0.52 


10* 


2 


230 


2 


700 


2.85 


8.71 


51 


-0 


550 


1 


840 


2.55 


2.40 


92 


-1 


600 


1.820 


1.80 


0.50 


11* 


1 


932 





425 


3.60 


0.40 


52 


-0 


600 


3 


020 


3.30 


2.29 


93 


-1 


580 


0.020 


1.80 


9.77 


12* 


1 


950 





750 


6.00 


21.82 


53 


-0 


650 


3 


150 


6.30 


2.75 


94 


-1 


550 


-0.080 


1.80 


28.18 


13* 


2 


150 


2 


650 


2.10 


2.51 


54 


-0 


637 


1 


820 


2.55 


16.60 


95 


-1 


720 


2.900 


3.00 


2.29 


14* 


2 


100 


2 


732 


2.25 


0.12 


55 


-0 


637 


1 


680 


4.50 


1.82 


96 


-1 


580 


-0.040 


7.95 


151.84 


15*t 


2 


020 


3 


332 


5.25 


3.78 


56 


-0 


750 


-2 


175 


4.35 


0.44 


97 


-1 


650 


0.040 


2.85 


10.71 


16* 


1 


350 





250 


2.55 


1.83 


57 


-0 


680 





070 


3.00 


2.19 


98 


-1 


620 


-0.050 


2.10 


4.68 


17* 


1 


312 





200 


2.55 


0.69 


58 


-0 


670 





562 


4.35 


34.67 


99 


-1 


620 


0.220 


2.85 


1.91 


18 


1 


250 


-1 


037 


5.55 


2.86 


59 


-0 


760 


-0 


450 


4.35 


1.32 


100 


-1 


720 


0.280 


6.00 


25.12 


19* 


1 


467 


1 


867 


5.55 


6.18 


60 


-0 


780 


1 


852 


4.35 


17.72 


101 


-1 


750 


-0.100 


3.60 


63.73 


20* 


1 


267 





267 


2.10 


0.36 


61 


-0 


780 





400 


2.55 


12.59 


102 


-1 


900 


1.870 


3.60 


1.15 


21* 


1 


380 


1 


460 


3.60 


0.79 


62 


-0 


800 


-1 


160 


4.80 


1.00 


103 


-1 


820 


-0.250 


4.20 


21.06 


22* 


1 


360 


1 


650 


3.30 


0.52 


63 


-0 


800 


1 


900 


2.70 


14.07 


104 


-1 


840 


0.000 


3.90 


23.19 


23 


1 


080 


-1 


560 


2.85 


0.25 


64 


-0 


880 





600 


6.75 


15.14 


105 


-1 


800 


-0.350 


2.85 


0.79 


24 


1 


100 


-1 


300 


6.90 


11.71 


65 


-1 


050 


2 


350 


2.25 


1.74 


106 


-1 


880 


-0.160 


14.40 


52.48 


25* 


1 


250 


1 


120 


3.30 


2.00 


66 


-1 


000 


-1 


620 


3.30 


2.05 


107 


-1 


820 


-1.800 


6.90 


1.20 


26 


1 


050 


-1 


150 


6.30 


17.72 


67 


-1 


050 





125 


3.30 


5.87 


108 


-1 


920 


-0.250 


3.00 


2.51 


27 


1 


000 





160 


3.30 


0.24 


68 


-1 


120 





560 


0.90 


0.03 


109 


-2 


000 


0.450 


3.30 


0.24 


28 





820 


-i 




5.10 


1.82 


69 


-1 
-i 


iUU 


-i 


oUU 


4.05 


7.43 


110 


1 
-i 


OOO 


-1.880 


3.60 


0.60 


29 





900 


-0 


667 


3.60 


0.36 


70 


-1 


160 


1 


950 


2.70 


2.00 


111 


-2 


020 


0.190 


3.90 


15.22 


30 





900 





050 


2.55 


0.17 


71 


-1 


160 





750 


Z.55 


1.51 


112 


-2 


180 


1.950 


Z.85 


n no 

0.93 


31 





800 





400 


4.80 


9.12 


72 


-1 


240 


2 


840 


4.05 


9.55 


113 


-2 


050 


0.320 


1.80 


0.26 


32 





900 


2 


080 


6.30 


3.31 


73 


-1 


160 





650 


2.55 


0.46 


114 


-2 


200 


1.425 


3.00 


1.58 


33 





667 





160 


3.30 


0.28 


74 


-1 


200 





150 


3.90 


5.35 


115 


-2 


400 


2.867 


6.75 


3.88 


34 





680 


1 


960 


4.50 


2.75 


75 


-1 


180 


1 


767 


2.10 


2.63 


116 


-2 


380 


2.020 


4.80 


7.24 


35 





440 





000 


4.95 


10.96 


76 


-1 


200 


1 


650 


5.25 


39.81 


117 


-2 


532 


-0.680 


3.60 


15.69 


36 





550 


1 


960 


2.10 


1.91 


77t 


-1 


420 


2 


350 


9.60 


37.84 


118 


-2 


600 


-0.732 


2.85 


2.49 


37 





562 


1 


900 


1.80 


1.51 


78 


-1 


300 


1 


780 


2.85 


3.80 


119 


-2 


780 


0.980 


2.10 


0.10 


38 





562 


1 


820 


2.70 


1.28 


79 


-1 


300 


1 


850 


3.90 


1.91 


120 


-2 


820 


1.040 


4.50 


0.55 


39 





400 


-0 


160 


3.90 


2.51 


80 


-1 


300 


-0 


532 


3.60 


0.38 


121 


-3 


000 


1.000 


7.95 


1.74 


40 





000 


-1 


820 


2.10 


0.13 


81 


-1 


350 


-0 


300 


4.95 


55.39 


122 


-3 


320 


0.960 


4.05 


0.55 


41 





050 





132 


6.75 


50.42 


82 


-1 


467 


1 


732 


3.30 


2.40 















^Lucke-Hodge number, see §2.1. 

'^OS'sets in degrees from 5*^ 24™, —69° 50' (B1950.0), increasing north and west. 

■^Radii in arcmin, sec §2.1. 

'^lO"*^ erg s~^ cm~^ 

^Irregular, nonspherical morphology. 

*Within observed WISP field. 
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Table A2 

Cavity parameters for Monte Carlo models. 



ShelP 


xb 




radius'^ 


LMC-r 


2.025 


3.325 


24.0 


LMC-2 


-2.225 


-0.475 


30.0 


LMC-3 


-1.075 


0.350 


34.5 


LMC-4 


-1.100 


2.375 


40.5 


LMC-5 


-0.600 


2.900 


27.0 


LMC-6* 


1.950 


0.575 


21.0 


LMC-7* 


2.125 


-0.250 


27.0 


LMC-8 


1.075 


-1.225 


30.0 


LMC-9 


-0.725 


-1.775 


30.0 



'"From Mcaburn (1980), sec §2.2. 

b Offsets in degrees from 5'' 24'", 
-69° 50' (B1950.0), increasing north 
and west. 

'^Radii in arcminutcs, sec §2.2. 
*Within observed WISP field. 



